Virial Theorem and Universality in a Unitary Fermi Gas 
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Unitary Fermi gases, where the scattering length is large compared to the interparticle spacing, 
can have universal properties, which are independent of the details of the interparticle interactions 
when the range of the scattering potential is negligible. We prepare an optically-trapped, unitary 
Fermi gas of 6 Li, tuned just above the center of a broad Feshbach resonance. In agreement with the 
universal hypothesis, we observe that this strongly-interacting many-body system obeys the virial 
theorem for an ideal gas over a wide range of temperatures. Based on this result, we suggest a 
simple volume thermometry method for unitary gases. We also show that the observed breathing 
mode frequency, which is close to the unitary hydrodynamic value over a wide range of temperature, 
is consistent with a universal hydrodynamic gas with nearly isentropic dynamics. 

PACS numbers: 03.75.Ss, 32.80.Pj 
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Universal behavior is believed to approximately de- 
scribe a variety of strongly- interacting Fermi systems, 
such as neutron stars Q, lllS, and resonantly inter- 
acting atomic Fermi gases Universal Fermi systems 
satisfy the unitary condition, where the zero energy scat- 
tering length greatly exceeds the interparticle spacing, 
while the range of the scattering potential is negligible. 
Unitary conditions are produced in an optically trapped 
Fermi gas by using a magnetic field to tune near 
a broad Feshbach resonance, where strong interactions 
are observed [6j- According to the universal hypothe- 
sis, the interparticle spacing sets the only natural length 
scale, and the system can exhibit universal thermody- 
namics If unitary Fermi gases satisfy 
this hypothesis, they can be used to test predictions in 
fields well outside of atomic physics. 

Universality at low temperature has been tested in 
measurements of two-component, strongly-interacting 
Fermi gases. The ratio, (3, of the interaction e nergy to the 
local Fermi energy has been measured [a, |lfj, |vjj |l2t |l3j 
and is believed to be a zero-temperature, universal 
parameter |lj. The measurements are in reasonable 
agreement with recent predictions [Til Hfij . Universal- 
ity also may explain the small line shifts observed in 
radio-frequency spectroscopy measurements Uni- 
versality at low temperature is further supported by 
the spatial profile 0, 0, 0] and breathing mode fre- 
quency 0,0,0 of a trapped, highly degenerate, uni- 
tary Fermi gas. However, for a Fermi gas tuned to a 
broad Feshbach resonance, there have been no model- 
independent tests of the universal hypothesis over a wide 
temperature range, and these tests impact all theoretical 
predictions of the thermodynamics. 

In this Letter, we show theoretically and experimen- 
tally that a harmonically trapped, unitary 6 Li Fermi gas 
obeys the virial theorem for an ideal gas over a wide 
range of temperature, as predicted by universal thermo- 
dynamics. Although the gas is a strongly-interacting, 
many-body system, we find that mean square width of 



the cloud varies linearly with the total energy, as pre- 
dicted by the virial theorem. We also show that uni- 
versal hydrodynamics, under isentropic conditions, may 
describe the behavior of the observed breathing mode fre- 
quency which is near the zero temperature unitary 
hydrodynamic value over a wide range of temperatures, 
and has no mean field shift. 

We test the universal hypothesis in a highly degener- 
ate, unitary Fermi gas of 6 Li, confined in a CO2 laser 
trap. This far-detuned trap provides a nearly harmonic 
potential, which is the same for all atoms, both paired 
and unpaired. We employ a 50-50 mixture of the two 
lowest spin-up and spin-down hyperfine states, and tune 
a bias magnetic field to 840 G, just above the center of a 
broad s-wave Feshbach resonance. Forced evaporation in 
the optical trap is then used to cool the gas, as described 
previously 0, 0, 000. After evaporation, we ob- 
tain a total of N = 2.0(0.2) x 10 5 atoms, and the cloud 
has a nearly zero-temperature Thomas-Fermi profile, as 
expected for a unitary gas at zero temperature 0|. 

From the measured trap frequencies we obtain tu±_ = 
t/lJJJ^ =2irx 1696(10) Hz, uj x /uj y = 1.107(0.004), and 
A = lo z /uj±_ = 0.045. The typical Fermi temperature 
(at the trap center for a noninteracting gas) is Tp — 
(SNy^h^xLUyUJz) 1 / 3 /k B ^ 2.4 (jK, small compared to 
the final trap depth of Uo/ks — 35 [iK. The coupling 
parameter of the strongly-interacting gas at B — 840 G 
is kpa ~ —30.0, where Tikp = \J1m ksTp is the Fermi 
momentum, and a — a(B) is the zero-energy scattering 
length estimated from the measurements of Bartenstein 
et al., HJ. 

We now show that universality requires such a 
strongly-interacting, unitary Fermi gas to obey the virial 
theorem for a harmonically-trapped ideal gas at all tem- 
peratures. According to the universal hypothesis, the 
thermodynamic functions which describe the gas must 
be independent of the interaction parameters and can de- 
pend only on the total density n and temperature T ||. 

Consider first the local energy AE (kinetic and inter- 
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action energy) contained in a small volume AV of gas 
centered at position x in a harmonic trap. Assume that 
the volume AV contains a fixed number of atoms AN, 
so that n = AN/AV, where J d 3 xn(x) = N is the total 
number of trapped atoms. 

For such a unitary gas, the local energy must be of the 
general form, 



AE = ANe F (n)f E 



T 



T F (n) 



(1) 



Here, the natural energy scale for atoms of mass m 
is taken to be kBTp(n) — e F (n) with e F (n) = 
ti 2 (3n 2 n) 2 / 3 /(2m). With this definition, e F (n) is the 
local Fermi energy corresponding to the density n and 
T F (n) is the corresponding local Fermi temperature. 
Note that for a zero-temperature ideal Fermi gas, we 
have /e = 3/5, for a zero- temperature unitary gas f E = 
3(1 + p)/5, while for a classical gas, f E = (3/2)T/T F (n). 
The corresponding local entropy AS takes the form 



AS = ANk B f s 



T 



T F {n) 



(2) 



where k E fs is the average entropy per particle, which 
can contain normal and superfluid contributions. 

The local pressure of the gas is readily determined from 
the relation P = -[d(AE)/d(AV)] A N,AS- From Eq. EJ 
we see that holding the local entropy constant requires 
fs = constant, which in turn means that we hold the lo- 
cal reduced temperature constant in taking the derivative 
of AE with respect to volume AV. Hence, we need only 
to find the volume derivative of the local Fermi energy, 
which yields the local pressure, 



P=-£(n,T), 



(3) 



where the local energy density (total kinetic plus 
interaction energy per unit volume) is £(n,T) — 
ne F (n) f F [T/T F (n)]. Eq. prelates the pressure and local 
energy density for the unitary gas in the same way as for 
an ideal, noninteracting homogeneous gas, although the 
energy densities are quite different. Eq. [3] for a unitary 
gas was obtained previously in Ref. 

In mechanical equilibrium, the balance of the forces 
arising from the pressure P and trapping potential U 
yields 



VP(x) + n(x) V£/(x) = 0. 



(4) 



Taking an inner product of Eq. 0] with x and using 
x • Vt/(x) = 2 U (x) for a harmonic trap, one readily ob- 
tains N(U) = (3/2) Jd 3 xP(x), where (U) is the aver- 
age potential energy per particle. Using f d 3 x£(x) = 
E - N(U) and Eq.^Jthen yields 
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FIG. 1: Verifying the virial theorem in a unitary Fermi gas of 
6 Li: (x 2 }/{x 2 (0)} versus E/Eo showing linear scaling. Here 
(a; 2 ) is the measured transverse mean square size. E is the 
total energy, calculated as in Ref. H3- E and (a; 2 (0)} denote 
ground state values. 



Hence, universality requires a unitary Fermi gas to obey 
the virial theorem for an ideal gas. Since the mean square 
size is oc (U) , Eq. [S] is equivalent to 



E 



(x 2 {E)) 
(x 2 (P )) E Q 



(6) 



which we use to verify the theorem. Here Eq is the ground 
state energy of the cloud. 

Energy is first added to gas, always starting from the 
lowest temperatures, by abruptly releasing the cloud and 
then recapturing it after a short expansion time theat E3 • 
During the expansion time, the total kinetic and inter- 
action energy is conserved. When the trapping potential 
U(x) is reinstated, the potential energy of the expanded 
gas is larger than that of the initially trapped gas, in- 
creasing the total energy to E (theat), which is a known 
function of theat • After waiting for the cloud to reach 
equilibrium, the sample is released from the trap. The 
mean square width (x 2 ) is estimated by fitting a one- 
dimensional, finite-temperature, Thomas-Fermi profile to 
the spatial distribution of the cloud 0] , which is imaged 
after a fixed expansion time of 1 ms. 

Fig. n shows (x 2 ) as a function E. The dashed line 
shows the fit, (x 2 )/(x 2 (0)) = 1.03 (0.02) E/E , which is 
in close agreement with the virial theorem prediction of 
Eq.Elfor a unitary gas. Eq.[5]is therefore verified because 
N(U(T = 0)) = E /2, which follows generally from the 
equation of state for a unitary gas at zero temperature 0, 
[l(J,[ljj. The Fermi radius a' x of the trapped unitary cloud 
is measured from a fit at nearly zero temperature 



(5) 

and determines (x 2 (T = 0)) = a /2 /8 as well as E /N 
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3mul(x 2 {T = 0)). 

From these results, we see that despite the strong, 
many-body interactions, the total potential energy of the 
unitary gas is half of the total energy. Hence, the sum of 
the kinetic and interaction energies must be exactly half 
of the total energy at all temperatures. 

We suggest that an empirical thermometry method can 
be based on the virial theorem result for a unitary Fermi 
gas. One can simply measure (x 2 ) for a unitary gas near 
a broad Feshbach resonance. This determines the total 
energy according to Eq. El i.e., since E/N = 3muj 2 (x 2 ). 
Universality then requires that E/Eq is in one-to-one cor- 
respondence with the reduced temperature, T/Tp. 

An approximate empirical reduced temperature can be 
determined by assuming that E/Eq obeys ideal gas scal- 
ing with reduced temperature. At a later time, the em- 
pirical temperature can be calibrated, either theoretically 
or experimentally. 

A theoretical calibration, relating the empirical and 
theoretical reduced temperatures, can be accomplished 
using an exact calculation of E(T/T F )/ Eq for a unitary 
gas. While this calibration method is not useful for mea- 
surements of the energy versus temperature (where it 
is a tautology), it can be used to estimate temperature 
for precise quantitative comparisons between predictions 
and measurements of condensed pair fractions 12 lL |22| , 
the gap [23|], collective mode damping rates etc. 

An experimental consistency check of the calibration 
can be done by measuring the total entropy of the uni- 
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tary gas versus energy, S(E). E is known from 



S 



can be determined, albeit in a model-dependent way at 
present Ha 1241 . using an adiabatic sweep of the mag- 
netic field |2^| from the unitary regime to a weakly in- 
teracting regime either well below I22L |23| or well above 
resonance |2_J , where S is calculated |l3j, |24| . Then, the 
fundamental relation, 1/T = (dS/dE) yields the tem- 
perature of the unitary gas for a given E, which can be 
compared to the empirical temperature measured at the 
corresponding E. 

We now consider the implications of universal hydrody- 
namics under isentropic expansion conditions, which may 
apply to the radial breathing mode of a unitary Fermi 
gas We find that the damping rate reveals a tran- 
sition in behavior, while the frequency remains close to 
the zero-temperature unitary hydrodynamic value over a 
wide range of temperature [ijj . 

Under locally isentropic conditions, the stream veloc- 
ities of the normal and superfluid components must be 
equal (since the entropy per particle is different for the 
superfluid and normal components). Then, we can as- 
sume that the total density n(x, t) and the stream veloc- 
ity u(x, t) obey a simple hydrodynamic equation of mo- 
tion. The convective derivative of the stream velocity u is 
the local acceleration which depends on the forces arising 
from the local pressure P and the trap potential U. For 
irrotational flow, V x u = 0, we have u • Vu = V(u 2 /2). 



du _ /m o \ \/f 



VP 



(7) 



where m is the atom mass |25| . 

For a unitary Fermi gas, the local pressure, Eq. Intakes 
the general form P(ra,T) = n 5 / 3 f P [T/T F (n)]. 

Initially, the gas is contained in a harmonic trap at a 
uniform temperature To and has a density no = no(x), 
where x is the position in the initial distribution. The 
initial pressure Po(x) — n^ 3 fp[T /Tp(n )]. Force bal- 
ance, Eq.0J requires VxPo(x)/n — — Vsf7(x). 

The hydrodynamic equation of motion can be solved by 
assuming a scaling ansatz |26l \'2l\ , where each dimension 
changes by a scale factor &j(t), i = x,y,z, and &i(0) = 
1. The density and stream velocity then take the forms 
n(x, t) = no(x)/r and u x — xb x (t)/b x (t) and similarly 
for u y ,u z . Here x = (x/b x ,y/b y ,z/b z ) is the position at 
time t — for an atom which is at position x at time t 
and r = b x b y b z is the volume scale factor. The scaling 
ansatz is exact if the gas is contained in a harmonic trap 
and the pressure takes the form P = en 1 , where c and 7 
are constants [2(3, |2?| ■ 

At nonzero temperature, the pressure does not in gen- 
eral obey such a simple power law, as the function f P can 
be dependent on T/Tp(n) in a complicated way. How- 
ever, for a gas expanding under isentropic conditions, 
an exact scaling solution can be obtained, and it pre- 
dicts temperature-independent expansion. The results 
correspond closely with our measurements, where nearly 
temperature-independent breathing frequencies are ob- 
served, as we now show. 

If the gas is locally isentropic, according to Eq. [2 the 
local reduced temperature does not change as the gas 
expands, i.e., T(x)/T F [n(x)] = T /Tp[no(x)]. Using n — 
uq/T then requires T(x) = T0/T 2 / 3 . If local equilibrium 
is maintained, the pressure P is then simply related to Pq. 
Using f P [T/T F [n(x)}] = fp[T /T F [n (5c)}] we must have 
P(x) = P (x)/T 5 / 3 . Then with x = b x x and V x -> V x in 
Eq.|lfor P (x), we obtain VP(x)/n(x) = -V[/(x)/r 2 / 3 . 
Using this result in Eq. [7] yields 



du 



171 2 TT, \ U(X) 

-u 2 + [/(x) - ^ 



(8) 



where U (x) is the trapping potential while U (x) arises 
from the pressure force, which is evaluated at x = x/b x (t) 
and similarly for y, z. Eq.|Hlcan be obtained for a stro ngly 
collisional gas by using a phase-space scaling ansatz |2j| • 
We see that under isentropic conditions, the gas obeys 
the same hydrodynamic equation at all temperatures as 
for a zero temperature gas, where P = cn 5 ' 3 and c is 
a constant, even though the pressure and density may 
assume a more complicated form. For a harmonic po- 
tential, all terms in Eq. |H1 are linear in x and the scaling 
solution is exact [2(|. Hence, after release from a har- 
monic trap, where £/(x) — > in Eq.|Hl the cloud expands 
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FIG. 2: Radial breathing frequency uj/u± versus empirical 
temperature f for a unitary gas of 6 Li from Ref. |T^I . showing 
small variation. Open circles: Measured frequencies; Solid 
dots: Data after correction for anharmonicity; Lower dot- 
dashed line: unitary hydrodynamic frequency, lu/uj± — 1.84. 
Upper dashed line: noninteracting gas frequency 2uj x /uj± — 
2.10. 



precisely by a scale transformation at all temperatures, 
similar to that observed in our previous measurements 
for a strongly-interacting Fermi gas at very low temper- 
ature 

The breathing frequencies are readily obtained using 
b t = 1 + £j, where e.; << 1. For the radial mode in a 
cylindrically-symmetrically trap, i.e., with u> a 
loj_ >> u) z , one obtains uj = ^/l0/3 wj_ = 
However, for our trap conditions, the exact result is 
lu = 1.84 independent of the mean field and super- 



— UJy — 

1.83 wj.. 



fluid contributions which are included in the general form 
of the energy density. 

Fig. [5] shows the measured breathing mode frequen- 
cies as a function of empirical temperature T, 0.11 < 
T < 1.54. Here, T is determined from the measured 
spatial profiles as in R ef. Il3l| and can be calibrated to 
the theoretical profiles [KJ, which show that the corre- 
sponding T/Tp varies from 0.12 to 1.1. After correction 
for anharmonicity, the frequency is close to the universal 
hydrodynamic value over the range of temperatures stud- 
ied. These results are consistent with nearly isentropic 
conditions at the highest temperatures, although the sys- 
tem is likely to be changing from a superfluid to a uni- 
tary collisional fluid, as suggested by the transition in the 
damping rate [lif. The observed hydrodynamic behav- 
ior at higher temperatures is not explained by existing 
theories. In particular, the momentum relaxation rate 
predicted in two-body collision models is much too small 
to explain the observed hydrodynamic behavior 29] . 

We are indebted to Jason Ho for a discussion of isen- 
tropic hydrodynamics, and to Jason Ho, Kathy Levin 
and Qijin Chen for a critical reading of the manuscript. 
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